PROBABILISTIC MODELING OF THE RENAL 
STONE FORMATION MODULE 


Integrated Medical Model 
Exploration Medicine Capabilities Project 


Lauren M. Best 1 

Jerry G. Myers 1 
Debra A. Goodenow 1 
Michael P. McRae 2 
Travis C. Jackson 3 


1 NASA Glenn Research Center, Cleveland, OH 

2 Rice University, Houston, TX 

3 University of Nebraska, Lincoln, NE 


www.nasa.gov i 


Outline 

• Overview of I MM 

• Renal Stone Risk Module 

• Probabilistic Wrapper - Front End 

• Biochemistry Model 

• Probabilistic Wrapper - Back End 




www.nasa.gov 2 



Integrated Medical Model (I MM) 



• IMM is a probabilistic tool - quantifies the probability and consequences of 
medical risks 


• Identifies medical resources necessary to optimize health and mission 


success 


Potential Medical 
Condition 



Evaluate 
with IMM 




Likelihood of 
occurrence , probable 
severity of occurrence, 
and optimization of 
treatment and resources 
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Renal Stone Risk Module 



• Astronaut urine samples have 
exhibited characteristics associated 
with renal stone formation, 
demonstrating an increased risk that 
clinically significant renal stones 
could form in flight 

• External module to IMM 

- Adds higher fidelity to model where little 
data is available 

• Deterministic Model ( Kassemi et. a!) 
- integrates physics and physiology 
into formation of renal stone based 
on probabilistic input 

• Probabilistic Wrapper - feeds 
needed information into deterministic 
model and transforms output into a 
useful distribution of risk of renal 
stones 
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Probabilistic Wrapper: Front-End 

Through Monte Carlo methods, randomly samples array of probability distributions 
to compile a set of random inputs into biochemistry and deterministic model 



Monte Carlo 
Simulation 
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Distributions of total 
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Monte Carlo Sampling 


Lognormal probability distributions 

formed from x and a and number of 
desired MC samples (N) in MATLAB 

- Truncated to the 95 th Cl 

- Example: 


Pea ~ 


% 

x t 


Ca-i 
Ca i' 


Parameters: 

• Calcium 

• Oxalate 

• Uric Acid 

• Citrate 

• Sulfate 

• Phosphate 



Magnesium 

Sodium 

Potassium 

Ammonium 

Creatinine 

(Creatine) 


Xca N . 

• Parameters randomly sampled and 
placed into matrix 

• Parameter matrix is divided by 
randomly sampled volume to create 
total concentrations (mol/L) 



[Pea P()x Pcreat 
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Example histogram of Calcium (truncated) 
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Biochemistry Model 



• JESS Biochemistry model - compiles 100+ applicable 
urine chemistry reactions and computes outputs based on 
system parameters and thermodynamics 

• Transforms total concentration matrix, M, via system of 
equations into free ion concentrations (Cj) based on urine 
chemistry reactions 

• Free ion must be converted to activity (a;) for deterministic 
model: 


di = , Yi = species activity coefficient 

C 

c e = standard concentration 
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Biochemistry Model 

• Relative Supersaturation calcium oxalate (RS Ca0x ) 

- More of the dissolved material than could be dissolved by 
the solvent 

- RS: Metric that represents the ability of spontaneous 
crystallization of the solution 

RS — [Q?aQ?xKcaOx] ^ 

• Necessary outputs to deterministic model: activity 
(indirect), pH, temperature, pressure, RS Ca0x 
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Probabilistic Wrapper: Back-End 

Couples deterministic model output and randomly sampled input 
parameters to quantify the risk of stone formation and treatment 

Utilizes two major sections: (Initially design by Mike McRae and Travis Jackson) 

1 . Probability of stone formation 

2. Probability that a stone passes spontaneously (without treatment) 




Probability distribution of male stone-forming astronauts 



Probability of Stone Formation 
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Back-End Continued... 

Stone Formation: 

- Determined by Multivariate logistic regression equation below via 
MATLAB statistical functions 

- Additional parameters can be added with additional 3 factors 

^ ^ 1 + e -(«+PilRS]+p 2 [^\+P 3 [TV]+p 4 {pH]+(] 5 lMg]) 

Stone Passage: 

- Determined via spontaneous passage rates related to renal stone 
width and ureter stone location 

- Treatment determined by the probability of stone passage: 

P(Treatment ) = 1 — P(Passage ) 
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